function y = boltz(x, p1, p2, p3)
y = p1./(1+exp(-(x-p2)*p3));
return;
